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Abstract 



The relative smoothing rates of various gauge field smoothing algorithms 
are investigated on 0(a 2 )-improved SU{2>) Yang-Mills gauge field configura- 
tions. In particular, an 0(a 2 )-improved version of APE smearing is motivated 
by considerations of smeared link projection and cooling. The extent to which 
the established benefits of improved cooling carry over to improved smearing 
is critically examined. We consider representative gauge field configurations 
generated with an 0(a 2 )-improved gauge field action on 16 3 x 32 lattices at 
(3 = 4.38 and 24 3 x 36 lattices at (5 = 5.00 having lattice spacings of 0.165(2) 
fm and 0.077(1) fm respectively. While the merits of improved algorithms are 
clearly displayed for the coarse lattice spacing, the fine lattice results put the 
various algorithms on a more equal footing and allow a quantitative calibra- 
tion of the smoothing rates for the various algorithms. We find the relative 
rate of variation in the action may be succinctly described in terms of simple 
calibration formulae which accurately describe the relative smoothness of the 
gauge field configurations at a microscopic level. 
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I. INTRODUCTION 



Gauge field smoothing algorithms are now widely used in lattice gauge theory studies as 
effective tools for constructing operators providing enhanced overlap between the vacuum 
and the hadronic state under investigation. APE smearing is now widely used in creating 
improved operators for static quark potential studies, or creating orbitally excited and hybrid 
mesons from the vacuum. Studies of perfect actions have lead to the construction of "Fat 
Link" fermion actions in which the links appearing in the fermion action are APE 

smeared. Such actions display better chiral behavior and reduced exceptional configuration 
problems. 

Both cooling and smearing algorithms have been used extensively in studies of QCD 
vacuum structure, where the lattice operators of interest suffer from large multiplicative 
renormalizations |9,l0]. Here the suppression of short distance physics is key to removing 
these perturbative renormalizations. 

Unimproved smoothing algorithms such as standard cooling |TT|-[T^| or standard APE 
smearing |l| introduce significant errors on each sweep through the lattice. These errors 
act to underestimate the action [HJ and spoil instantons as the action falls below the one- 
instanton bound. The problems may be circumvented by adding additional irrelevant oper- 
ators to the action tuned to stabilize instantons [^,|T^]. 

Such improved cooling algorithms are central to studies of topology and instantons in the 
QCD vacuum. There thousands of sweeps over the lattice are required to evolve a typical 
gauge field configuration to the self-dual limit. It is well established that the use of improved 
algorithms is central to achieving the required level of accuracy. 

In this paper we introduce an 0(a 2 )-improved form of APE smearing and examine the 
extent to which the benefits of improvement in cooling algorithms carry over to improved 
smearing. To carefully examine this new algorithm we create gauge field configurations with 
an (9(a 2 )-improved gauge action. For this investigation we select Symanzik improvement for 
the gauge action |T7|]. We consider two sets of gauge field configurations; a coarse 16 3 x 32 
lattice at (3 = 4.38 with a ~ 0.165(2) fm, and a fine 24 3 x 36 lattice at f3 = 5.00 providing 
a ~ 0.077(1) fm. 

These lattices are sufficiently fine that we expect similar results for other choices of ac- 
tion improvement schemes such as Iwasaki []TSJ| or DBW2 [[TJJ, actions explored in Ref. ||20 



We also define an 0(a 2 )-improved non-abelian field strength tensor and construct the cor- 
responding improved topological charge operator. 

While the merits of improved algorithms are clearly displayed for the coarse lattice 
spacing, the fine lattice results put the various algorithms on a more equal footing. Moreover, 
on the fine lattice we no longer witness transitions between topological charge values as a 
function of smoothing sweeps. 

Finally, we calibrate the relative smoothing rates of standard cooling, APE smearing, 
improved cooling and improved smearing using the action as a measure of the smoothness. 
The action is selected as it varies rapidly under cooling. The action evolution has only a mild 
configuration dependence H] allowing the consideration of only a few configurations in the 
calibration process. We focus our computational resources on numerous cooling schemes, 
including an improved version of APE smearing referred to as improved smearing. In all we 
consider fourteen smoothing algorithms for two hundred sweeps on seventeen configurations. 
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This is computationally equivalent to the more standard study of one or two algorithms on 
the order of a hundred configurations. This calibration analysis is also an extension of an 
earlier analysis []55| which focused on unimproved algorithms. 



The plan of this paper is as follows: Section || describes the lattice action used in this 
simulation. The improved field strength tensor and associated topological charge operator 
are described in Section |HI[ In motivating improved smearing we begin in Section [IV| with 
a brief review of improved cooling followed by our improved smearing algorithm. Section [V] 
presents the results of our numerical simulations. The calibration results are discussed in 
Section [VT] and a summary of the findings is given in Section |V1J| , where the connection to 
previous studies drawing relations between cooling and physical properties |9|,^|^]] is made. 

II. LATTICE GAUGE ACTION 

The tree-level O (a 2 ) -improved action is defined as, 

Sg = f £ ^Tr(l - P^ix)) - -§- 2 £ ReTr(l - Rfax% (1) 

U>Ll V>Ll 

where P^ u and R^ u are defined as 

Pfax) = Ufa) Ufa + A) Ul(x + Z)Ufa), (2) 
Rfax) = Upix) U v {x + fjt) U v (x + z> + /2) Ufa + 20) Ufa + v) Ufa) 

+ Ufa) Ufa + p) Ufa + 2/2) Ufa + ft + 0) Ufa + u) Ufa). (3) 



The link product Rfax) denotes the rectangular 1x2 and 2x1 plaquettes. u is the tadpole 
improvement factor that largely corrects for the quantum renormalization of the coefficient 
for the rectangles relative to the plaquette. We employ the plaquette measure for the mean 
link 

n \ 1/4 

u = i^-ReTr (Pfax))j , (4) 

where the angular brackets indicate averaging over x and fi ^ v. 

Gauge configurations are generated using the Cabibbo-Marinari [25| pseudo-heat-bath 



algorithm with three diagonal SU(2) subgroups cycled twice. Simulations are performed 
using a parallel algorithm on a Thinking Machines Corporation CM-5 with appropriate link 
partitioning . 

Configurations are generated on a 16 3 x 32 lattice at (5 — 4.38 and a 24 3 x 36 lattice 
at /3 — 5.00. Configurations are selected after 5000 thermalization sweeps from a cold 
start, and every 500 sweeps thereafter with a fixed mean-link value. Lattice parameters are 
summarized in Table ffl. 
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III. TOPOLOGICAL CHARGE OPERATOR 



The topological charge of a gauge field configuration provides a particularly sensitive 
indicator of the performance of various smoothing algorithms. The topological charge is 
related to the field strength tensor by 



9 2 

Q = Y1 l( x ) = zl ^r^W^ Tr (F pu {x)F pa {x)) , 



(5) 



where q(x) is the topological charge density. An expression for may be obtained by 
expanding the definition of the Wilson loop. Consider a loop C in the [i-v plane 



Cfj,u(x) = "Pexp (ig j> A(x) ■ dx^j 

l + ig(£ A(x) - dx )~^ (f A ( x ) ■ dx ) + °(i 



(6) 



The line integral is easily evaluated using Stokes theorem and a Taylor expansion of d fl A u (x) 
about xq 



A(x) ■ dx 



F/j,v 

+\ (xlDl + xlD 2 u ) i^(x ) + 0(a 2 g 



2 a 4 ) 



(7) 



The integration limits are determined by the size of the Wilson loop. Positioning the ex- 
pansion point xq at the center of the Wilson loop one finds 



j£ xi A{x) ■ dx = a 2 F^(x ) + ^ {pi + D 2 ) F^(x ) + ... 
I A(x) ■ dx = 2a 2 F flu (x ) + (4Dj + D 2 ) F^(x ) + ... 

J 2 x 1 X _ 

<f A(x) ■ dx = 2a 2 F llu (x ) + ^(Dl+ AD 2 \ F^{x ) + .... 



(8) 
(9) 
(10) 



Hence, F pu may be extracted from consideration of the lxl plaquette alone. To isolate 
the second term of the expansion in Eq. (P), one takes advantage of the Hermitian nature 
of F^x). Constructing F pu (x) symmetrically about x leads to 



gF pu (x) = — 



0${x) - 0<M{x)) - -Tr (p${x) - 0% 



)t, 



x 



(11) 



TABLE I. Parameters of the numerical simulations. 



Action 


Volume 


A^Thcrm 


A^Samp 





a (fm) 


u 


Physical Volume fm 


Improved 


16 3 x 32 


5000 


500 


4.38 


0.165(2) 


0.8761 


2.64 3 x 5.28 


Improved 


24 3 x 36 


5000 


500 


5.00 


0.077(1) 


0.9029 


1.848 3 x 2.772 
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(1) I u J x l 



FIG. 1. Graphical representation of the link products summed in creating 0$(x). 



where O^J(x) is the sum of 1 x 1 Wilson loops illustrated in Fig. [T], 

0$(x) = Ufa) Ufa + A) Ufa + v) Ul(x) 

+ Ufa) U^(x + v — p) Ul(x — p) Ufa — p) 

+ Ufa - p) Ufa - p-u) Ufa -p-u) Ufa - v) 

+ Ufa - v) Ufa - v) Ufa + p-i>) Ufa) . (12) 

This well known definition for Ffax) is commonly employed in the Clover term of the 
Sheikholeslami-Wohlert |27j improved quark action. Although, not everyone enforces the 
traceless nature of the Gell-Man matrices by subtracting off the trace as in Eq. ([TT|) . 

Unfortunately, this definition has large 0(a 2 ) errors. These errors are most apparent in 
the topological charge. Even after hundreds of sweeps of cooling this simple definition fails 
to take on integer values. Errors are typically at the 10% level. 

To improve the topological charge operator, we improve the definition of the field strength 
tensor F^ v by removing 0(a 2 ) errors with a linear combination of plaquette and rectangle 
Wilson loops 

0$(x) = <*0${x) + %lj$(x), (13) 

Uq 

where the tadpole coefficient u is defined in Eq. (£|). Here lj?J(x) is the link products of 
1x2 and 2x1 rectangles in the [i-v plane 

Xf}{x) = Ufa)Ufa + jl)Ufa + 2p)Ufa + p + P)Ufa + u)Ufa) 
+ Ufa)Ufa + p)Ufa + p + P)Ufa + 2P)Ufa + i>)Ufa) 
+ Ufa)Ufa + P)Ufa - p + 2P)Ufa - p + v)Ufa - jl)Ufa - p) 
+ Ufa)Ufa -p + v)Ufa -2p + P)Ufa - 2p)Ufa - 2p)Ufa - p) 
+ U*fa — p)UUx — 2p)U\{x — 2/2 — P)Ufa — 2p — P)Ufa — p — P)Ufa — 0) 
+ Wfa — p)Ul(x — p — u)Ul(x — p — 2P)U fa — p — 2P)Ufa — 2v)Ufa — P) 
+ Ufa - P)Ufa - P)Ufa + p - i>)Ufa + 2p- P)Ufa + p)Ufa) 
+ Ufa - P)Ufa - 2v)Ufa - 2v)Ufa + P- 2v)Ufa + p- P)Ufa) , (14) 

depicted in Fig. The coefficients c\ and C2 are determined by combining Eqs. (||), @ 
and (|10D to remove 0(a 2 ) errors. 
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FIG. 2. Graphical representation of the link products summed in creating O f _ a )(x). 



lxl 



A(x) ■ dx 



f A(x) ■ dx + f A(x) ■ dx 

J2xl Jlx2 



indicating the coefficients for the 0(a 2 )-improved F^ v of Eq. ([13]) are 



5 A 1 

ci = - and c 2 = --. 

6 O 



(15) 



(16) 



IV. COOLING AND SMEARING ALGORITHMS 

In this section we motivate and introduce an improved version of APE smearing. The 
motivation is based on improved cooling, and therefore we begin with a very brief overview 
of cooling algorithms. Cooling consists of minimizing the local action effectively one link at 
the time. The local action is that contribution to the action associated with a single link. 

A. Improved Cooling 

Standard cooling minimizes the local Wilson action associated with a link at each link 
update. The local action associated with the link U^(x) is proportional to 



Si{x,n) = ReTr(l - U ll (x)E IJ3/ (x)) 



(17) 



V 



where S MJy (x) is the sum of the two staples associated with U^x) lying in the [i-v plane 

= U v {x + p)Ul(x + v)Ut{x) + Ul(x + /t - u)Ul(x - v)U v {x - v) . (18) 

The local action associated with the link U^(x) is minimized in standard cooling by a mini- 
mization process, i.e., by replacing the original link by the link U^{x) which optimizes 



maxReTr U^x) ^ 



x 



(19) 



V 
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Improved cooling proceeds in exactly the same manner, but with the plaquette-based 
staples replaced with the linear combination of plaquette-based staples and rectangle-based 
staples. For improved cooling we use 



max Re Tr 



V 



pu V^) i 



(20) 



V 



Where 



yj - 5 v , 1 yR 
^ 3 12u 2 



(21) 



and 



= U v {x + (a; + v + £)£/*(£ + 2j>)[/J( a; + 0)U}(x) 
+ Upix + fi)U u {x + 2fr)Ul{x + A + i>)t/J(a; + *>)£#(a0 
+ £7,(x + p)Ul(x + i>)Ul(x + z> - fi)Ul(x - /t)Z7 M (x - //) 

(a; — z>) 

+ f/ M (x + p)Ul{x + 2/t - z>)l^(z + /} - u)Ul(x - P)U v (x - v) 
+ ?7^(x + /t — fyUUx — v)Jj\{x — v — jl)U v {x — v — p)UJx — ft) . 



(22) 



The mean-field factor uq is updated following each sweep through the lattice and rapidly 
goes to 1 as perturbative tadpole contributions are removed. As we wish to study 0(a 2 )- 
improved cooling, the coefficients are based on those of the action |jPT| , pS|| . While this action 
improves contact with QCD, this choice of action does not completely stabilize instantons 
on the lattice [|T5|,[^J. Removal of 0(a 2 ) errors and stabilization of instantons requires the 
consideration of additional loops |T6| , |30| . The preferred algorithm for finding the U^(x) which 
maximizes Eq. ( |19D is based on the Cabbibo-Marinari |25] pseudo-heat-bath algorithm for 
constructing 577(3) -color gauge configurations. There, operations are performed at the 
SU(2) level where the algorithm is transparent. 

An element of SU(2) may be parameterized as, U = a^I + ia ■ a, where a is real and 
a 2 = 1. Since sums of products of SU(2) matrices are proportional to SU(2) matrices, 



where U^x) G SU(2) and 



E KM 

V 



k 2 = det 



(23) 




(24) 



The maximum of the expression 



ReTr ^^E^ 



x 



v 



ReTr (kU^x) U^x) 



(25) 
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is achieved when 




Re Tr (U^x)!! „(x) j = Re Tr(J) , (26) 
which requires the link to be updated as 



At the 577(3) level, we successively apply this algorithm to the three diagonal 577(2) sub- 
groups of 577(3) 

When considering an improved cooling algorithm it is crucial to loop over sufficient SU(2) 
subgroups to ensure that the subtle effects of the higher dimension operators introduced in 
improving the action are reflected in the final SU(3) link. It is easy to imagine that only two 
or three SU(2) subgroups may not take the SU(3) link close enough to the optimal link for 
the effects of improvement to be properly seen. We find consideration of the three diagonal 
SU(2) subgroups looped over twice to be optimal. In fact we have seen round off errors 
actually increase the action if too many loops are made. 

The Mo factor of Eq. (EHI) is not held fixed during the improved cooling iteration. Starting 
from the value determined during the thermalization process, uq is updated after every sweep 
through the lattice. After a few sweeps of improved cooling the value quickly converges to 
1 which is a good indication that short distance tadpole effects are being removed in the 
smoothing procedure. 

Finally it must be understood that cooling proceeds effectively one link at a time. That 
is, links involved in constructing T, 1 (x) are not updated simultaneously with U^(x). In fact 
it is a nontrivial task identifying which links can be simultaneously updated on a massively 
parallel computer p6| . 



B. Improved Smearing 

In this section we consider the APE smearing [|TJ algorithm. We extend this algorithm 
to produce an improved version, motivated by the success of the improved cooling program. 



1. Reunitarization of the Links 

APE smearing is a gauge equivariantQ |32| prescription for smearing a link U^x) with its 
nearest neighbors U^{x + z>), where v is transverse to (x. The APE smearing process takes 
the form 



Gauge equivariance means that if two starting gauge configurations are related by a gauge 
transformation then the respective smeared configurations are also related by the same gauge 
transformation. 
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U,(x) — %{x) = (1 - a)U,(x) + | ? St^x) , (28) 

= P^(x) , (29) 

where T,^ v {x) is the sum of the two staples associated with U^(x) defined in Eq. and V 
is a projection operator which projects the smeared link back onto SU(3). 

The reunitarization procedure is of central importance when smearing is applied to the 
gauge links because the APE smearing operation produces links outside the SU(3) gauge 
group in the intermediate stage of smearing. The link between smearing and cooling is 
established via the manner in which smeared links are projected back onto SU(3). The 
preferred approach is to select the link variable U^{x) that maximizes the quantity 



Re Tr(^ (*)[/;(*) J • (30) 
In this case a clear connection to cooling is established when a — > 1 as 

maxReTr (u^(x)U'^(x)^ =maxReTr U^{x) ^2 ^^(x) , (31) 



which is precisely the condition of Eq. (19) for cooling. In other words, projection of the 
smeared link back onto the gauge group via Eq. ([H]) selects the link which minimizes the 
local action. Improvement of the staple based on the action will aid in removing 0(a 2 ) 
errors encountered in the SU(3) projection. In practice, the same Cabibbo-Marinari-based 
cooling method of operating on SU(2) subgroups may be used to obtain the ultimate SU(3) 
link. 

However there is one very important difference remaining between APE smearing and 
cooling. While cooling effectively updates one link at a time, feeding the smoothed link 
immediately into the next link update, APE smearing proceeds uniformly with all links 
being simultaneously updated. Smoothed links are not introduced into the algorithm until 
the next iteration of the APE smearing process takes place. 

It is this latter point that provides the form factor interpretation of APE smearing in 
fat-link fermion actions M. There smearing can be understood to introduce a form factor 
suppressing the coupling of gluons to quarks at the edge of the Brillouin zone where lattice 
artifacts are most problematic. The form factor analysis restricts the smearing fraction to 
the range < a < 3/4. Indeed in practice, smearing fractions beyond 3/4 do not lead to 
smooth gauge configurations. 



2. Improving Smearing 

Having made firm contact with cooling it is clear that replacing the simple staple of 
Eq. (|Tj^) with the improved staple of Eq. fl2T| ) may lead to a realization of the benefits of 
improved cooling within APE smearing. Hence we define improved smearing to be an APE 
smearing step in which 
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U,(x) — Ul(x) = (1 - a)U,(x) + | £ EjJW 



(32) 
(33) 



where S^(a; 



is the improved staple of Eq. (|21|) . 
Signatures of improvement include the preservation of structures in the action density 
distribution under smearing. 0(a 2 ) errors in the standard Wilson action act to underesti- 
mate the local action and destroy topologically nontrivial field configurations ]T3] . Improved 
smearing will have reduced 0(a 2 ) errors and hence better preserve topologically nontrivial 
field configurations. Hence an associated signature of improvement is the stability of topo- 
logical charge under hundreds of smearing sweeps. Of course, one could alter the coefficients 
of the improvement terms to stabilize instantons |T5| at the expense of introducing 0(a 2 ) 
errors into the smoothing action. 

The extended nature of the staple will alter the stability range of a. Where standard 
APE smearing provides the range < a < 3/4, we find the upper stability limit lies below 
a = 0.6. 



V. NUMERICAL SIMULATIONS 

We analyze two sets of gauge field configurations generated using the Cabibbo- 
Marinari pseudo-heat-bath algorithm with three diagonal SU(2) subgroups looped over 
twice. Details may be found in Table |I|. As discussed in the introduction, analysis of a few 
configurations proves to be sufficient to resolve the nature of the algorithms under investiga- 
tion. We consider eleven 16 3 x 32 configurations and six 24 3 x 36 configurations. For each 
configuration we separately perform 200 sweeps of cooling and 200 sweeps of improved cool- 
ing. We explore 200 sweeps of APE smearing at seven values of the smearing fraction and 
200 sweeps of improved smearing at five values of the smearing fraction. For APE smearing 
we consider a = 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, and 0.70. Similarly for improved smearing 
we consider a = 0.10 to 0.50 at intervals of 0.10. The extended nature of the staple alters 
the stability range of a to lie below a = 0.6. 

For clarity, we define the number of times an algorithm is applied to the entire lattice 
as n c , riic, n ape (a) and ni ape (a) for cooling, improved cooling, APE smearing and improved 
smearing respectively. We monitor both the total action normalized to the single instanton 
action So = 8n 2 /g 2 and the topological charge operators, Ql and Q J ™ P , from which we 
observe their evolution as a function of the appropriate sweep variable and smearing fraction 
a. 



A. The Influence Of The Number Of Subgroups On The Gauge Group. 

In this section we describe the influence of including additional SU(2) subgroups in con- 
structing the gauge group SU(3) and explore the impact it has on the smoothing procedure. 
The Cabibbo-Marinari algorithm []25| constructs the SU(N) gauge groups using SU(2) sub- 
groups. It is understood that the minimal set required to construct SU(3) matrices is two 
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180 200 



FIG. 3. The evolution of the topological charge estimated by the improved operator as a func- 
tion of standard cooling sweeps n c for various numbers of SU{2) subgroups. The curves are for 
a typical configuration from the 16 3 x 32 lattices where a = 0.165(2) fm. The parameter cycle 
describes the number of times the three diagonal SU (2) subgroups are cycled over. 



diagonal 577(2) subgroups. After having performed cooling on gauge field configurations 
we noticed that the resulting cooled gauge field configurations were not smooth even after 
a large number of smoothing steps. Adding a third SU{2) subgroup made a significant 
difference in the resulting smoothness. 

We explored further by simply performing additional cycles, n cyc i e , around the three 
diagonal 577(2) subgroups. We monitored the smoothing rate using both the Standard 
Wilson and the improved action according to the cycle number being set to n cyc i e = 1,2 and 
3. Based on the evolution of the action, we found that the optimum cooling rate on gauge 
field configurations is achieved using the three diagonal 577(2) subgroups cycled over twice, 
^cycie = 2. Cycling more than twice provides very little further reduction of the action and 
round off errors may actually increase the action on occasion. Hence two cycles over the three 
diagonal SU(2) subgroups is sufficient to precisely create the SU(3) link which minimizes 
the local action. This determination is crucial to ensuring the effects of our improved action 
are fully reflected in the 577(3) link. 

We also monitored the evolution of the topological charge with respect to the above 
number of cycles. On the 16 3 x 32 lattice with spacing of a = 0.165(2) fm, we observe a 
disagreement of the trajectories for the topological charge for different numbers of the 577(2) 
subgroups cycles. Fig. |3] displays results for standard cooling and Fig. |] displays similar 
results for improved cooling. A comparison of Figs. |3| and |] indicates the trajectories also 
differ between cooling and improved cooling. 

Hence we see subtle differences in the algorithms leading to dramatic differences in the 
topological charge. One must conclude that a lattice spacing of 0.165(2) fm is too coarse for 
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20 40 60 80 100 120 140 160 180 200 

Ic 

FIG. 4. The evolution of the topological charge estimated by the improved operator as a func- 
tion of improved cooling sweeps n\ c for various numbers of cycles over the three diagonal SU(2) 
subgroups. The curves are for the same configuration illustrated in Fig. ^. 

a serious study of topology in SU(3) gauge fields using the algorithms considered here. The 
difficulty lies in the fact that the algorithms have a dislocation threshold typically a little over 
two lattice spacings |TB| , PT| . Instantons with a size smaller than the dislocation threshold are 
removed during the process of smoothing. While this threshold has the desireable property of 
removing lattice artifacts, we note that twice the lattice spacing is 0.33 fm. Minor differences 
in the dislocation thresholds of the various algorithms will cause some (anti) instantons to 
survive under one algorithm where they are removed under another. 

On the other hand, a lattice spacing the order of 0.077(1) fm appears to allow a mean- 
ingful study of topology in SU(3) gauge theory. Twice this lattice spacing is 0.15 fm, well 
below the typical size of instantons [ f2~4"| . 

We also note here the accuracy with which our improved topological charge operator 
reproduces integer values. These results should be contrasted with the usual 10% errors of 
the unimproved operator at similar lattice spacings. Such errors on a topological charge of 
5 can lead to the uncomfortable result of Q ~ 4.5 when the unimproved operator is used. 

With the finer 24 3 x 36 lattice at f3 = 5.00, we observe perfect agreement among trajecto- 
ries for different numbers of cycles of the three SU(2) subgroups. Moreover, the topological 
charge remains stable for hundreds of sweeps following the first three sweeps. Figures |5| 
and H] compare the topological charge evolution for cooling versus improved cooling for six 
configurations. In every case, the two algorithms produce the same topological charge for a 
given configuration. 
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FIG. 5. The evolution curve for the topological charge estimated via the improved operator as 
a function of cooling sweeps n c for six configurations on the 24 3 x 36 lattices at = 5.00 where 
a = 0.077(1) fm. 
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FIG. 6. The evolution curve for the topological charge estimated via the improved operator as 
a function of improved cooling sweeps n\ c for the same six configurations from the 24 3 x 36 lattices 
at P = 5.00 illustrated in Fig. |. 
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FIG. 7. The ratio S/Sq as a function of standard cooling sweeps n c for five configurations on 
the 24 3 x 36 lattice at f3 = 5.0. The single instanton action is Sq = 8ir 2 /g 2 . 

B. The Action 

We begin by considering the action evolution on both lattices. Here we report the action 
divided by the single instanton action Sq = 8ir 2 /g 2 . It is important to note that although 
the 24 3 x 36 lattice has almost four times more lattice sites than the 16 3 x 32 lattice, the 
physical volume is smaller by almost a factor of three. As such the typical topological charges 
encountered are smaller in magnitude. 

Figs. ^, ||, |§ and ^ report the typical evolution of the action under standard cooling, 
improved cooling, APE smearing and improved smearing respectively. Inspection of the 
figures reveals that improved cooling preserves the action better than standard cooling over 
a couple hundred sweeps. As expected, standard APE smearing remains slower than cooling 
or improved cooling even at our most efficient smearing fraction (a = 0.70). Similar results 
are observed for improved smearing at our most efficient smearing fraction of a = 0.50 in 



Fig. 10 



Based on these observations, one concludes that the fastest way to remove the short range 
quantum fluctuations on an 0(a 2 ) gauge field configuration, is through standard cooling, 
which lowers the action more rapidly than improved cooling as a function of cooling sweep. 
In turn we see that improved cooling is faster than the maximum stable standard APE 
smearing, which is faster than the maximum stable improved smearing. This is illustrated by 
Figs. |?Hl0[ It is important to emphasize that the fastest way of removing these fluctuations 
is not necessarily the best as far as the topology is concerned. It is already established that 
the 0(a 2 ) errors of the standard Wilson action act to underestimate the action [ 15| . These 
errors spoil instantons which might otherwise survive under improved cooling. 
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FIG. 8. The ratio S/Sq as a function of improved cooling sweeps n\ c for five configurations on 
the 24 3 x 36 lattice at (3 = 5.0. The rate of cooling is seen to be somewhat slower than that for 
the standard cooling. 
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FIG. 9. The ratio S/Sq as a function of APE smearing sweeps n SLpc (a) for one configuration 
on the 24 3 x 36 lattice at [3 = 5.0. Each curve has an associated smearing fraction a. The rate of 
lowering the action for the maximum stable smearing fraction (~ 0.75) is seen to be less than that 
for the other standard or improved cooling. 
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FIG. 10. The ratio S/Sq as a function of improved smearing sweeps ni ape (a) for one configu- 
ration on the 24 3 x 36 lattice at /3 = 5.0. Each curve has an associated smearing fraction a. We 
see that this is the slowest of the four algorithms for lowering the action as a function of the sweep 
number. 



C. Topological Charge from Cooling and Smearing 



We begin by considering the 16 3 x 32 lattices having a lattice spacing of a = 0.165(2) 
fm. In Fig. [XT], we plot the evolution curve for the improved topological charge as a function 
of the cooling sweep number, n c , for six of our configurations. Similarly in Fig. 0, for the 
same six configurations we plot the improved topological charge operator but this time as 
a function of improved cooling sweep n\ c . The line types in Figs. |TT] and [12] correspond to 
the same underlying configurations and are to be directly compared. For example, the solid 
curve corresponds to the same gauge field configuration in both figures but with a different 
algorithm applied to it. From these two figures we notice the two cooling methods lead to 
completely different values for the topological charge. 

Improved cooling brings stability to the evolution of the topological charge whereas 
standard cooling gives rise to numerous fluctuations to the topological charge. For improved 
cooling, plateaus appear after about forty sweeps and persist for hundreds of sweeps. This 
is a celebrated feature of improved cooling. 

This algorithmic sensitivity of the topological charge is also seen in Figs. [13] and |T^, for 
APE and improved smearing on a single configuration (the solid line of Figs. |TT] and |12|) . 
Within APE smearing or improved smearing, the topological charge trajectories follow simi- 
lar patterns but different rates for various smearing fractions. However, APE smearing leads 
to values for the topological charge which are different from that obtained under improved 
smearing. An important point is that improved smearing stabilizes the topological charge 
at 95 sweeps for a = 0.5, whereas standard smearing shows no sign of stability until 140 
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FIG. 11. Q L mp versus n c for six configurations on the 16 3 x 32 lattices at (3 
fm. Each line corresponds to a different configuration. 
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versus n\ c for six configurations on the 16 x 32 lattices. The different line types 



identifying different configurations match the configurations identified in Fig. [11. 
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FIG. 13. The evolution of Qi^ p using APE smearing as a function of APE smearing sweep 
(a) on the 16 3 x 32 lattice at [5 = 4.38. Here different line types correspond to different 



80 100 120 140 160 180 200 
n (a) 

ape v ' 



n 



ape 



smearing fractions. 



iterations at a = 0.5. Hence we see significant improvement in the topological aspects of 
the gauge field configurations under improved smearing. 

On our finer lattice, we find a completely different behavior for the topological charge 
evolution. The topological charge is established very quickly; after a few sweeps in the case 
of cooling or improved cooling as illustrated in Figs. [5] and || The topological charge persists 
without fluctuation for hundreds of sweeps, both for cooling and for smearing as illustrated 
in Figs. [15| and [16| for APE and improved smearing respectively. Moreover, the topological 
charge is independent of the smearing algorithm. 

These results are to be compared with Fig. |ll| and Fig. [12], where transitions are observed 
even with improved cooling on the coarser lattice with a = 0.165(2) fm. The results on our 
fine lattice suggest the characteristic size of instantons is much larger than the dislocation 
threshold, such that the topological structure of the gauge fields is smooth at the scale of 
the threshold. 

In Fig. [T5| for standard APE smearing we observe a slower convergence to integer topo- 



logical charge than in Fig. 16 for improved smearing when 0.1 < a < 0.5. This feature 



of improved smearing is illustrated in detail in Fig. |17|. However, APE smearing has the 
advantage to allow values for the smearing fraction up to a = 0.70 which cannot be accessed 
by improved smearing. 

Having demonstrated that it is possible to precisely match the behavior of the algorithms 
on fine lattice spacings, we proceed to calibrate the efficiency of these algorithms in the 
following section. 
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FIG. 14. The evolution of Q^ ip using improved smearing as a function of APE smearing sweep 
^iape(oj) on the 16 3 x 32 lattice at f3 = 4.38. Here different line types correspond to different 
smearing fractions. 
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FIG. 15. The evolution of Q^ np using APE smearing as a function of APE smearing sweep 
^ape(a) on the 24 3 x 36 lattice at (3 = 5.00. Here different line types correspond to different 
smearing fractions. 
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FIG. 16. The evolution of Q^ ip using APE smearing as a function of APE smearing sweep 
^iape(«) on the 24 3 x 36 lattice at (5 = 5.00. Here different line types correspond to different 
smearing fractions. 
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FIG. 17. The evolution of the improved topological charge, Q^ p , as a function of standard 
APE smearing sweeps, n apc (a), for 0.1 < a < 0.7 (solid lines) is compared to improved smearing 
sweeps, ni a pe(a), (dotted-dashed lines) for the same smearing fractions 0.1 < a < 0.5 on the 
24 3 x 36 lattice at (3 = 5.00. The horizontal dotted-dashed line is Q l ™ p = -4. 
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VI. SMOOTHING ALGORITHM CALIBRATION 



Here we calibrate the relative rate at which quantum fluctuations are removed from 
typical field configurations by the various algorithms. The calibration is done using the 
action normalized to the single instanton action, S/Sq, on both, the 16 3 x 32 lattices and 
24 3 x 36 lattices. The action normalized to the single instantons action, S/ So, is of particular 
interest because it provides insight into the lattice content as well as the rate at which the 
quantum fluctuations are removed. 

While there is no doubt that the algorithms may be accurately calibrated on the fine 
24 3 x 36 lattice, the 16 3 x 32 lattice with a = 0.165(2) fm presents more of a challenge. As 
a result, in most cases we will show the graphs produced from the 16 3 x 32 lattice analysis 
and simply present the numerical results for both the 16 3 x 32 and 24 3 x 36 lattices. 

The numerical results are summarized in Tables [0], |TTTj, and |VI] for the 16 3 x 32 lattices 
and in Tables |TV|, [V], and [VTT| for the 24 3 x 36 lattices. 



A. APE Smearing and Improved Smearing Calibration. 

To calibrate the rate at which the algorithm reduces the action we record the nearest 
number of sweeps required to reach a given threshold in S/Sq. The action thresholds are 
spaced logarithmically to obtain a uniform distribution in the number of sweeps required to 
reach a threshold. The relative rates of smoothing are established by comparing the relative 
number of sweeps required to reach a particular threshold. 

Here we calibrate the APE smearing algorithm characterized by the smearing fraction 
a and the number of smearing iterations n apc (a). The different threshold crossings are 
characterized by the number of sweeps required to reach that threshold, n ape (a). In Fig. [L8| 
we show the number of sweeps required to reach a threshold when a = 0.5, n ape (0.50), 
relative to that required for other a values, n ape (a). We plot these relative smoothing 
rates as a function of n ape (a) such that low S/So thresholds are reached after hundreds of 
iterations of the smoothing algorithm. Fig. [19] shows similar results for improved smearing. 
In these figures and in the following analysis, we omit thresholds that result in fewer than 
five smoothing iterations as these points produce integer discretization errors of more than 
20%. 

Both standard and improved smearing algorithms have a relative smoothing rate which 
is independent of the amount of smoothing done. By calculating the average value for each 
of the bands in Figs. |T^ and [H^, we can investigate the dependence of the average relative 
smoothing rate (n ape (0.50)/n ape (a)) on a. 

Fig. ^0] illustrates a linear fit to the data constrained to pass through the origin. We find 
(n ape (0.50)/n ape (a)) = 2.00 a such that 

in agreement with our earlier analysis ||22j| . The extent to which this relationship holds can 
be verified |22| by plotting the ratio a'n&p^a') / an apc (a) and comparing the results to 1. 
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16 3 x 32 lattice at (3 = 4.38. From top to bottom the data point bands correspond to a 
0.5, 0.4, 0.3, 0.2, and 0.1. 
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FIG. 19. The ratio ni ape (0.50)/ni ape (a) versus ni ape (a) for numerous S/So thresholds on the 
16 3 x 32 lattice at (3 = 4.38. From top to bottom the data point bands correspond to a = 0.5, 0.4, 
0.3, 0.2, and 0.1. 
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FIG. 20. Illustration of the dependence of (n ape (0.50)/n a p e (a)) for APE smearing on the smear- 
ing fraction a. The solid line is a linear fit to the data constrained to pass through the origin. 



In plotting the band averages for the improved smearing algorithm of Fig. |19] in Fig. |21| , 
one finds a small deviation of the points from the line y = 2 a. This suggests that Eq. ( |34|) 
is not sufficiently general for the improved smearing case. 

A better approximation to establish the a dependence, that is similar to Eq. (0) and 
contains Eq. (K4T) is 



rciape(oQ _ ( Oi X " 
^lape(a) 



(35) 



where 5 is equal to one in the case of standard APE smearing and can deviate away from 
one for improved smearing. 

In Fig. |22|, the logarithm of Eq. (|3~5"D is plotted. The slope of the data provides 5 = 
0.914(1) for both the 16 3 x 32 and the 24 3 x 36 lattices. We also verified 5 = 1.00 for the 
APE smearing data. Fig. ^ plots the ratio of the left- and right-hand sides of Eq. fl35D as 
a function of ni ape (o;) for a' > a. The ratio is one as expected with 5% for large amounts 
of smearing where integer discretization errors are minimized. Throughout the following 
analysis, 6 is fixed at 0.914(1). 



B. APE and Improved smearing cross calibration. 



In this section we focus on the cross calibration of the smearing algorithms. In Fig. |24|, 
we compare the number of improved smearing sweeps required to reach a threshold for var- 
ious improved smearing fractions relative to APE smearing at a = 0.5. The lowest band 
corresponds to an improved smearing fraction of 0.10. From this, we conclude that for low 



23 




0.0 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

a/0.5 

FIG. 21. Illustration of the dependence of < rai ap c(0.50)/wi ap e(G!) > for APE smearing on the 
improved smearing fraction a. The solid line fit is constrained to pass through the origin. 
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FIG. 22. Illustration of the dependence of In (< ni ap e(0.50)/ni ape (a) >) on the improved smear- 
ing fraction a for improved smearing. The solid line fit indicates S = 0.914. 
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FIG. 23. Illustration of the degree to which the relation Eq. ( j35[) is satisfied for improved 
smearing. Here the entire data set is plotted for a and a' = 0.5, 0.4, 0.3, 0.2 and 0.1. Data are 
from 16 3 x 32 lattice at /3 = 4.38. 



a values APE smearing and improved smearing produce roughly equivalent smeared config- 
urations. However, there are some evident differences in the rate at which both algorithms 
perform. For intermediate to large a there is curvature in the bands. Early in the smear- 
ing process, fewer sweeps of improved smearing are required to reach a threshold. That is, 
improved smearing removes action faster than APE smearing in the early stages of smear- 
ing. This behavior is also manifest in the analogous results for the fine 24 3 x 36 lattice. 
As emphasized in the discussion surrounding Fig. [T7|, improved smearing also provides a 
topological charge closer to an integer than APE smearing. Together, these two properties 
of improved smearing identify a genuine improvement in the smearing process. 

For the coarse 16 3 x 32 lattice data, the bands are thick for large smearing fractions indi- 
cating improved smearing does perform significantly different from standard APE smearing. 
Contributions from individual configurations are clearly visible as lines within the bands. 
This structure is due to the coarse lattice spacing of 0.165(2) fm which reveals differences 
between the algorithms. Such structure is not seen in the fine 24 3 x 36 lattice results. There 
a precise calibration is possible. 

In Tables [[J and III, we report the averages of each band for APE and improved smearing 



on the 16 3 x 32 lattices. In Tables [TV| and |V| we report similar results for the 24 3 x 36 lattices. 
Based on equations Eq. ([34] ) and Eq. (|35|) for APE and improved smearing, we expect 



a n 



ape 



constant . 



(36) 



This ratio is plotted in Fig. ^5] where a rather mild dependence on ni ape (a) is revealed. 
Averaging these results provides 0.81(2) for the constant of Eq. (|3T)|). Similar results are 
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FIG. 24. The ratio n ape (0.50)/ni apc (a) versus ni apc (a) for numerous threshold actions on the 
16 3 x 32 lattice at (3 = 4.38. From top to bottom the data point bands correspond to improved 
smearing fractions a = 0.5, 0.4, 0.3, 0.2, and 0.1. 

TABLE II. The averages of the ratios < n a p C (0.50)/n ape (a) > and < ni ap e(0.50)/n ape (a) > for 
various smearing fractions a from the 16 3 x 32 lattice at (3 = 4.38. 
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a for APE smearing 
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nia pe (0.50) 0.227(1) 0.465(1) 0.706(1) 0.948(1) 1.189(1) 1.431(1) 1.673(1) 



TABLE III. The averages of the ratios < n apc (0.50)/ni ape (a) > and < ni ape (0.50)/ni apc (a) > 
for various smearing fractions a from the 16 3 x 32 lattice at (3 = 4.38. 
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TABLE IV. The averages of the ratios < n ap e(0.50)/n ape (a) > and < ni ape (0.50)/n ape (a) > 
for various smearing fractions a from the 24 3 x 36 lattice at (3 = 5.00. 
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TABLE V. The averages of the ratios < n a pe(0.50)/ni a pe(a) > and < ni a pe(0.50)/ni a pe(a) > 
for various smearing fractions a from the 24 3 x 36 lattice at (3 = 5.00. 









a for improved smearing 








0.10 


0.20 


0.30 


0.40 


0.50 


n apc (0.50) 


0.196(1) 


0.378(1) 


0.546(1) 


0.704(1) 


0.851(1) 


r*ia P c(0.50) 


0.228(1) 


0.442(1) 


0.641(1) 


0.826(1) 


1.0 



2.0 

l.B 

^ 1.6 
a 

1 1.4 

a 

/ 1.2 

S °- b 

/0.6 
b 0.4 
0.2 
0.0 

20 40 60 80 100 120 140 160 180 200 

n T (a) 

lape v ' 

FIG. 25. Illustration of the degree to which the relation Eq. ( ]36| ) is satisfied for calibration of 
the action under APE and improved smearing. Here the entire data set is plotted. 
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seen for the finer 24 3 x 36 lattice, but with greater precision in the calibration reflected in a 
narrower band. There the constant is also 0.81(2). However, it should also be noted that for 
a < 0.5, improved smearing achieves integer topological charge faster than standard APE 
smearing. 



C. Calibration of Cooling and Smearing 



In this section we apply the ansatz of equations Eq. ( p4|) and Eq. ( p5|) to relate the 
cooling and smearing algorithms. Fig. ^ displays results comparing cooling and standard 
APE smearing. For a as small as 0.1 it takes about 75 sweeps of APE smearing compared 
with 5 sweeps of cooling to arrive at an equivalent action. On the other end of the smearing 
fraction spectrum, we note the bands become very thick. 

The calibration of these ratios indicates 

U % , = 0.59(1), (37) 
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FIG. 26. The ratio ra c /n ape (a) versus n ape (a 
lattice at j3 = 4.38. From top down the data point bands correspond to a 
0.2, and 0.1. 



0.7, 0.6, 0.5, 0.4, 0.3, 



in agreement with that obtained in |22| where the analysis was performed on unimproved 
gauge configurations. The reduction in 0(a 2 ) errors in the gauge field action affect both 
algorithms similarly such that the calibration of the relative smoothing rates remains unal- 
tered. 

Further broadening of the bands is observed when comparing improved cooling with 



APE smearing as illustrated in Fig. The precision of improved cooling relative to APE 
smearing leads to very different smoothed configurations at this coarse lattice spacing of 
0.165(2) fm. This indicates the algorithms are sufficiently different, that an accurate and 
meaningful calibration is impossible. 



This effect is not observed when we pass to our fine lattice spacing as displayed in Fig. |28. 
We remind the reader that the thickness of the band for small numbers of smoothing sweeps 
is simply due to the ratio of small integers taken in plotting the y-axis values. 

The real test of improved smearing is the extent to which the algorithm can preserve 
action associated with topological objects and thus maintain better agreement with more 
precise algorithms including cooling and improved cooling. Fig. |2^ displays results for the 
calibration of improved cooling with improved smearing. Comparing these results for each 
smearing fraction, a, with that for improved cooling and standard smearing in Fig. ^7 
reveals that the improved smearing algorithm, which was seen to be better than standard 
APE smearing algorithm does not perform as well as the improved cooling algorithm. 

Similar results are seen in Fig. [3D] where standard cooling is compared with improved 
smearing. Hence the annealing of the links in the process of cooling, where cooled links are 
immediately passed into the determination of the next cooled link, is key to the precision 
with which cooling can preserve topological structure. 
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FIG. 28. 
lattice at (3 
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The ratio ni c /n ape (a) versus n ape (a) for numerous action thresholds on the 24 x 36 
= 5.00. From top down the data point bands correspond to a = 0.7, 0.6, 0.5, 0.4, 0.3, 



29 



0.5 




20 40 60 80 100 120 140 160 180 200 

n T (a) 

Iape v ' 

FIG. 29. The ratio ni c /ni ape (a) versus ni ape (a) for numerous action thresholds on the 16 3 x 32 
lattice at f3 = 4.38. From top down the data point bands correspond to a = 0.5, 0.4, 0.3, 0.2, and 
0.1. 
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FIG. 30. The ratio n c /ni apc (a) versus ni ape (a) for numerous action thresholds on the 16 3 x 32 
lattice at (5 = 4.38. From top down the data point bands correspond to a = 0.5, 0.4, 0.3, 0.2, and 
0.1. 
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TABLE VI. Calibration coefficients for various smoothing algorithms on the 16 3 x 32 lattice at 
(3 = 4.38. Entries describe the relative smoothing rate for the algorithm ratio formed by selecting 
an entry from the numerator column and dividing it by the heading of the denominator columns. 



For example equation Eq. (37) corresponds to the first column of the third row. 



Numerator 


an ape (a) 


Denominator 
a 6 riiape(a) 


n c 


nic 


a' n ape (a') 


1.00(2) 


0.81(2) 


1.69(3) 


1.30(2) 


a /<5 "-iape(a') 


1.25(3) 


1.01(2) 


2.13(5) 


1.61(3) 


n c 


0.59(1) 


0.47(1) 


1 


0.75(1) 


nic 


0.77(1) 


0.62(1) 


1.33(2) 


1 



TABLE VII. Calibration coefficients for various smoothing algorithms on the 24 3 x 36 lattice 
at P = 5.00. Entries describe the relative smoothing rate for the algorithm ratio formed by selecting 
an entry from the numerator column and dividing it by the heading of the denominator columns. 



For example equation Eq. (36) corresponds to the second column of the first row. 



Numerator 


a^apc(a) 


Denominator 
a 5 ni ap e(a) 


n c 


ni c 


a' n ap e(a') 


1.00(2) 


0.81(2) 


1.64(2) 


1.37(1) 


a /5 nia pe (a') 


1.25(3) 


1.01(2) 


2.04(4) 


1.67(3) 


n c 


0.611(9) 


0.49(1) 


1 


0.84(1) 


nic 


0.734(8) 


0.60(1) 


1.19(1) 


1 



Calibration of the smoothing rates as measured by the action for the algorithms un- 
der investigation are summarized in Tables [VI] and |V11| . The entries describe the relative 
smoothing rate for the algorithm ratio formed by selecting an entry from the numerator col- 
umn and dividing it by the heading of the denominator columns. The entry comparing APE 
smearing with itself reports the level to which the ansatz of equation Eq. (|34"D is satisfied. 
Similarly the entry comparing improved smearing with itself reports the level to which the 
ansatz of equation Eq. (|35|) is satisfied. 



D. Cooling versus Improved cooling. 



Figure |31] reports a comparison of standard cooling with improved cooling on eleven 
configurations from the coarse 16 3 x 32 lattice. There the ratio n c /rii c < 1 confirms the 
expectation that standard cooling does not preserve action on the lattice as well as the 
(9(a 2 )-improved cooling. Fewer standard cooling sweeps are required to reach the same 
action threshold. Calibration of the algorithms appears plausible for the first 80 sweeps of 
improved cooling, after which the two algorithms smooth the configurations in very different 
manners. Any calibration at this lattice spacing is only very approximate beyond 80 sweeps 
of improved cooling where distinct configuration-dependent trajectories become visible. This 
result is contrasted by the analogous analysis on our fine 24 3 x 36 lattice illustrated in 
Fig. While n c /ni c remains less than one, it is closer to one here than for the coarser 
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FIG. 31. The ratio n c /n\ c versus rvy c for numerous action thresholds on the 16 3 x 32 lat- 
tice at (3 = 4.38. The significant differences between the algorithms are revealed by the 
gauge-configuration dependence of the trajectories. 



lattice as one might expect. 



VII. CONCLUSIONS 

We have introduced an improved version of the APE smearing algorithm founded on 
the connection between cooling Eq. (|l~9"D and the projection of the APE smeared link back 
to the SU(3) gauge group via Eq. QHTj) ■ This connection motivates the use of additional 
extended paths combined with the standard "staple" as governed by the action to reduce 
the introduction of 0(a 2 ) errors in the smearing projection process. 

Clear signs of improvement are observed. For a given smearing fraction a defined in 
equation Q2"8|), improved smearing preserves the action better than standard APE smearing at 
each smearing sweep. At the same time improved smearing brings the improved topological 
charge to an integer value faster than standard APE smearing. 

The extended nature of the "staple" in improved smearing reduces the stability regime for 
the smearing fraction. We found the improved smearing algorithm to be stable for a < 0.5. 
At a = 0.6 the algorithm is unstable whereas standard APE smearing remains stable for 
a < 0.75. 

Given the wide variety of smoothing algorithms under investigation in the field of lattice 
gauge theory, we have cross calibrated the speed with which the algorithms remove action 
from the field configurations. In particular we have cross calibrated the smoothing rates of 
APE smearing at seven values of the smearing fraction; improved smearing at five values of 
the smearing fraction; cooling; and improved cooling. We explored smearing fractions in 0.1 
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FIG. 32. The ratio n c /n\ c versus n\ c for numerous action thresholds on the 24 3 x 36 lattice at 
P = 5.00. 



intervals starting at a = 0.1. 

The calibration has been investigated over a range of 200 sweeps for each smearing 
algorithm on 0(a 2 )-improved gauge field configurations. The results of this analysis allows 
one to make qualitative comparisons between cooling and smearing algorithms and in fact 
make quantitative comparisons of smearing algorithms with different smearing fractions on 
lattices as coarse as 0.165(2) fm. On our fine lattice where the lattice spacing is 0.077(1) 
fm, the calibration is quantitative in general. 

We have found the relative smoothing rates are described via simple relationships as 
reported in Tables VI and VII for our coarse 16 3 x 32 and fine 24 3 x 36 lattices respectively. 
There the sensitivity of the calibration results on the lattice spacing may be reviewed. A 
noteworthy point is that we discovered a necessary correction to the APE smearing ratio 



rule when improved smearing is considered. These algorithms may be calibrated via 



n ape {a ) _ a , /ii m [a ) / a 



nupe(a') 



^ape(a) ^lape(a) 



a' 



6 



and pcv ' = — (38) 



for APE smearing and improved smearing respectively. We find S = 0.914(1) without a 
significant dependence on the lattice spacing. 

Having cross calibrated these smoothing algorithms, we now proceed to make contact 
with physical phenomena |4],^,^4[] . In particular, we note that it is possible to build in a 



length scale beyond which cooling does not affect the links ||23|| . It would be interesting to 
explore these techniques in the context of APE and Improved Smearing. Using a random 
walk argument, one can postulate a cooling radius 

r CO o\ = Cy/n^a, (39) 
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where a is the lattice spacing and c is a constant independent of (3 f24" |. It has been shown 



that phenomena taken from simulation results with invariant a^fn~ c scale very well |[24| . The 
effective range for APE smearing has been estimated using analytic methods [|J. For small 
smearing fraction, a, the effective range is 

1 



r ape = -j=^an ape (a)a. (40) 

The product of a and n ape (a) defines r ape in agreement with the results presented here. 
Eq. ( |3lf ) indicates that this relation holds even for large a. Results of our analysis contained 
in tables |VI| and [VI I| allow one to link Eqs. (|39| ) and (|40| ) and thus determine the constant 



c. For sufficiently fine lattices c is argued to be independent of (3 MM and this is already 



supported to some extent by the similarity of the entries in Tables |\T]and [VI I[ For example, 
from Table |V1J[ , n c = 0.611(9) an ape («) such that 



^cooi = , =Jn~ c a = 0.739(5)^^0. 

'3(0.611(9)) 



The effective range for other smoothing algorithms may be derived from Eq. ([40l) in a 
similarly straight forward manner. 

Unfortunately a rigorous analysis of the scaling of the results of Tables [VI] and [VI I| is 
not possible. We have clear evidence that the topology of Yang-Mills gauge fields cannot be 
reliably studied using the algorithms presented here on lattice spacings as coarse as 0.165(2) 
fm. Different algorithms lead to different topological charges, differing quite widely in some 
cases as reported in Figs. [11] and 0. Moreover, subtle differences in the cooling algorithms 
can lead to different topological charge determinations as illustrated in Figs. [5] and |j. As 
discussed in Sec. V A, the proximity of the dislocation thresholds of the algorithms to the 



typical size of instantons and variations in the threshold from one algorithm to another 
causes some (anti) instantons to survive under one algorithm, whereas they are removed 
under another. 

In contrast, the fine 24 3 x 36 lattice results, where a = 0.077(1) fm, display excellent 
agreement among every smoothing algorithm considered. In this case it appears that the 
dislocation thresholds are smaller than the characteristic size of topological fluctuations^ 
such that the gauge fields are already sufficiently smooth to unambiguously extract the 
topology of the gauge fields. 

As a final comparison of the smoothing algorithms, we provide a visual representation of 



a gauge field configuration after applying various smoothing algorithms. Figure [33] illustrates 
a rendering of the topological charge density for a slice of one of the fine 24 3 x 36 lattice 
configurations. While our calibration has been carried out by considering the total action 
of the gauge fields, the following analysis allows us to examine the extent to which the 
calibration is accurate at a microscopic level. 

In Fig. red shading indicates large positive topological charge density with decreasing 
density becoming yellow in color, while blue shading indicates large in magnitude, negative 



2 We define a "topological fluctuation" to refer to objects with Q = ±1 but S/Sq > 1. 
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topological charge density decreasing in magnitude through the color green. Here cooling 
(a), improved cooling (b), APE smearing at a = 0.70 (c), APE smearing at a = 0.30 (d), 
improved smearing at a = 0.50 (e) and improved smearing at a = 0.30 (f), are compared 
at the number of smoothing iterations required for each algorithm to produce an approxi- 
mately equivalent smoothed gauge field configuration. While Fig. b) for improved cooling 
differs somewhat due to round off in the sweep number, the remaining plots compare very 
favorably with each other. These visualizations confirm that the different smoothing algo- 
rithms considered in this investigation can be accurately related via the calibration analysis 
presented here and summarized in Tables [VT] and |V11| . 
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FIG. 33. The topological charge density of a 24 3 x 36 lattice for fixed x coordinate. The 
instantons (anti-instantons) are colored red to yellow (blue to green). Fig. a) shows the topological 
charge density after 9 cooling sweeps. Each of the following figures display the result of a different 
smoothing algorithm calibrated according to Table |VII| to reproduce as closely as possible the 
results depicted in Fig. a). Fig. b) illustrates the topological charge density after 11 sweeps of 
improved cooling. Fig. c) shows the topological charge density after 21 APE smearing steps 
at a = 0.70. Fig. d) illustrates the topological charge density after 49 APE smearing steps at 
a = 0.30. In Fig. e) the topological charge density is displayed after 35 sweeps of improved 
smearing at a = 0.50. Finally, Fig. f) shows the topological charge density after 55 sweeps of 
improved smearing at a = 0.30. Apart from Fig. b) for improved cooling, which differs largely due 
to round off in the sweep number, all the plots compare very favorably with each other. 
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